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We investigate tlie accuracy of slow-roll inflation in light of current observational constraints, 
which do not allow for a large deviation from scale invariance. We investigate the applicability of 
the first and second order slow-roll approximations for inflationary models, including those with 
large running of the scalar spectral index. We compare the full numerical solutions with those given 
by the first and second order slow-roll formulae. We find that even first order slow-roll is generally 
accurate; the largest deviations arise in models with large running where the error in the power 
spectrum can be at the level of 1-2%. Most of this error comes from inaccuracy in the calculation 
of the slope and not of the running or higher order terms. Second order slow-roll does not improve 
the accuracy over first order. We also argue that in the basis eo — 1/H, e„+i = din |e„ |/dA'^, 
introduced by Schwarz et al. (2001), slow-roll does not require all of the parameters to be small. 
For example, even a divergent ea leads to finite solutions which are accurately described by a slow- 
roll approximation. Finally, we argue that power spectrum parametrization recently introduced by 
Abazaijan, Kadota and Stewart does not work for models where spectral index changes from red to 
blue, while the usual Taylor expansion remains a good approximation. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Inflation is a theory which postulates that a rapid ex- 
pansion of the universe occurred right after the Big Bang 
mialalj- Most inflationary models can be represented 
by an effective single field model with effective potential 
V. The inflaton with mass m rolls down the potential un- 
til the kinetic energy of the inflaton is greater than half of 
its potential energy. At this point the inflationary expan- 
sion of the universe stops and the next phase of reheat- 
ing occurs. During the inflationary expansion, the initial 
quantum fluctuations exponentially increase and become 
classical |E S H B IS Q • These classical fluctuations 
also seed the subsequent growth of large scale structure. 
There is a well defined procedure which allows us to find 
the spectrum of the fluctuations given the inflationary 
potential. Because exact solutions are numerically in- 
tensive several appoximations have been developed. The 
most common approximation is the slow-roll approxima- 
tion. Rece ntly the so-called uniform approximation was 
suggested 0[ 0. Reference developed improved 
WKB-type approximation. 

If the kinetic energy of the inflaton is much smaller 
than its potential energy, we say that the inflaton is 
slowly rolling down its potential. In this slow-roll ap- 
proximation we can obtain analytical formulae for the 
produced power spectrum in the form of a Taylor series 
expansion in a set of slow-roll parameters. The coeffi- 
cients in the Taylor expansion of the logarithm of the 
power spectrum in In k effectively define the slope — 1, 
running as and higher derivatives. We usually derive 
the slow-roll formulae through the time delay formalism 



or Bessel function approximation. Therefore there are 
some implied conditions on the accuracy of the slow-roll 
approximation depending upon the slow-roll parameters. 
References [l^, [13 found that there are areas in the slow- 
roll parameter space where the accuracy of slow-roll ap- 
proximation is questionable. This usually requires a large 
deviation of n., from 1. However it contradicts the latest 
observations [iIIitLIi^. 

Recently there has been a lot of renewed interest in 
models with large running of the scalar index [T^l20ll2]| . 
It is not clear whether slow-roll approximation is accurate 
in this area of parameter space, as in some expansions one 
of the slow roll parameters becomes large and the expan- 
sion is no longer well controlled Another issue is 
the question of where to stop the expansion. Although it 
is often assumed that the running is O ((ti^ — 1)^), ref- 
erence jl^ found that there are cases where it can be 
as large as rig — 1. In this case one should also consider 
the effect of including the running of the running of Us, 
i.e. the second derivative of Us over In A:. These are the 
issues addressed in this paper. We begin with a short 
review of the basic physics of inflation and the algorithm 
of numerical solutions to the inflationary equations, with 
more details given in appendix. We continue by compar- 
ing the numerical solutions to those given by slow-roll 
approximations and finally we present our conclusions. 

In this paper we use a standard convention for reduced 
Planck mass mpi — Gj^^'^ . 



II. INFLATIONARY BASICS 
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In the "Hamilton- Jacobi" formulation, the evolution of 
the Hubble parameter during inflation with poten- 
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tial V{(f>) is given by (e.g. see [23) 
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The number of e-folds N since some initial time is re- 
lated to the value of the scalar field (j) by 



dN _ 47r H{(l)) 



(2) 



'pi 



We will consider the situation when the value of the scalar 
field is growing in time, d(j)/dt > 0. Then by our conven- 
tion dN/dt is also positive, dN/dt > 0. 

In the literature, different sets of slow-roll parameters 
are used. Reference |23| introduce potential slow-roll pa- 
rameters which are constructed on the basis of the deriva- 
tives of the inflationary potential V{4>). Authors of 
define Hubble slow-roll parameters through the deriva- 
tives of the Hubble parameter H{(f)) with respect to the 
field 4> during inflation 
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In this parameterization when the inequality eH{4') < 1 
fails, the inflation immediately stops. Sometimes is 
also denoted £,h or though it can take negative values. 
In this paper we will use = ^^h ■ 

Reference [2^ introduces another basis of "horizon- 
flow" slow-roll parameters through the logarithmic 
derivative of the Hubble distance eo = dn = 1/H(N) 
with respect to the number of e-folds N to the end of 
inflation 
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The connection between any two of these sets can be 
found in e.g. |2^. Thus the first three horizon-flow slow- 
roll parameters are connected to the first three Hubble 
slow-roll parameters as 0, |2^ 



ei = en, 
£263 = Aeff -6eHi]H + 2^H- 
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There is an analytical connection between Hubble slow- 
roll parameters and potential slow-roll parameters |25| . 

References ^8, 29, 30^ 31j use differently defined sets 
of slow-roll parameters, but they still can be converted 
to the ones we have described here (e.g. see p6l|'). 

Thus any inflationary model can be completely de- 
scribed by the evolution of one of the sets of the pa- 
rameters. 



The condition for the inflation to occur is ei = e^f < 1 
or < 1 since ch = ey to first order. 

To find the power spectrum of the perturbations pro- 
duced by a single field inflation, one can follow the pre- 
scription of Grivell and Liddle ,3^] . One solves the equa- 
tion [Mill Hi 
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for each mode with wavenumber k and initial condition 
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-00. Then the spectrum of 



curvature perturbations is given by 



(11) 



The quantity z in equation (|10|l is defined as z = a(f)/ H 
for scalar modes and z — a for tensor modes. Then for 
scalar modes ^3 



-^ = 2a2i/2 
z dT'^ 
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One can parametrize the power spectrum of the scalar 
and tensor modes of the fluctuations amplified by the 
infiation as 
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around some conventional pivot point fc,. Leach et al. 

give expressions for the scalar spectral index , the 
running of the scalar spectral index a^, the tensor spec- 
tral index rit and the running of the tensor spectral index 
at in terms of the horizon-flow parameters. Here we re- 
produce their second order formulae for Ug — 1 and 

Us-l = -2ei - e2 - 2el 

-(2C7 + 3)eie2-Ce2e3, (14) 
as = -2eie2 - £263, (15) 

where C = jji + \n2-2^ -0.7296. 

Reference 15| also analyzes the accuracy of the approx- 
imation H13|) for parameterizing the inflationary power 
spectrum of fluctuations with /3 = for different values 
of the parameters r, Us and 

In this paper we will also compare the second order 
formulae H14I I15|) to the first order formulae given by 



Us - 1 = -2ei - 62, 
as = -2eie2 - 6263- 



(16) 
(17) 



The expression for as is the same as in the second order 
formula because the expression for as is derived using 
only first order expression for Ug. Thus the main differ- 
ence between first and second order formulae comes from 
the extra terms in the expression for — 1. 



3 



0.1 




-0.15 I ' ' ' ' ' 1 

0.8 0.9 1 1.1 1.2 

"s 

FIG. 1: The current constraints (68% and 95% confidence 
level contours) in the ns-a, plane from WMAP+SDSSgal 
(biKKer re d) and WMAP+SDSSlya (smaller, green) data jTa, 
Ha . l3a . l3q| . The constrained region clearly allows the value of 
the spectral scalar index Us to be around 1 and the running 
Os of the scalar spectral index to be significantly non-zero for 
either combination of the experiments. 



Current observational constraints on r, Ug and are 
given by [Hill 13- At 95% confidence level, the ten- 
sor to scalar ratio is i? < 0.50, which implies that the 
first horizon-flow parameter ei is much smaller than one. 
Current constraints on the scalar spectral index give us 
Ug = 0.98 ± 0.02, which in turn means that the second 
horizon-flow slow-roll parameter is much smaller than 
one. 

Present data does not require the presence of running 
in the primordial power spectrum |37i|. b ut running as 
large as ±0.03 is still allowed at 3-cr [ig- Regular in- 
flationary models usually predict \as\ ~ (n^ — 1)^ and 
so the running is of the order of 10"'^, as is the case for 
the minimally-coupled V{(j)) — Xcf)^ model with 60 e-folds 
remaining. 

But it is possible that |as| {ug — 1)^ and as < 0, 
which means that the main part in the running of the 
spectral index (|15|l is determined not by the first term 
— 2eie2, but by the second term —£2 £3- It happens when 
ksl ^ kil, and therefore there might be a situation when 

k3l>i. 

To summarize, if ris « 1 and ag is a small negative 
number, at some scale we might have £1 1, £2 ^ 1 
and I £3 1 > 1. Leach et al. |0| define inflation satisfying 
slow- roll under the condition |£„| <C 1, for all n > 0. In 
our case £3 > 1 , so the question arises as to whether slow- 
roll in this case is accurate or whether the approximation 
breaks down and one must also include terms with higher 
powers in £3. Does it mean that the inflation is not slow- 
roll and one must use full numerical solutions instead? 
And does it mean that one must also include the running 
of the running? These are the main questions we address 
in this paper. To address them we have developed the 
numerical code described in appendix A. 

How natural is it for inflation with a given number from 



50 to 70 e-folds remaining to produce a power spectrum 
with a changing tilt? In the absence of theoretical guid- 
ance on the inflationary space we cannot address this 
question simply. Authors of have produced about 
200,000 simulations of the inflationary flow equations for 
more or less "random" potentials, and calculated the ob- 
servable parameters {rig, ag, r, rit, at) of the resulting 
power spectra about 40 to 70 e-folds before the end of 
the inflation for each potential. About 80,000 of them fall 
into the area plotted on Fig. ^ Only the fifteen marked 
with larger yellow circles give a significant change in the 
tilt from red to blue, i.e. Ug ~ 1, a^ < —0.02. 

Choosing the Hubble parameter to be represented by a 
Taylor expansion in 4> with uniformly distributed coeffi- 
cients, as done in TP], does not necessarily correspond to 
the real infiationary priors [s^- We do not address this 
issue here; instead we want to simply stress that possi- 
bility of constructing a potential with a large running in 
the scalar power spectrum 40-70 e-folds before the end of 
the inflation exists. 



III. QUADRATIC POTENTIAL 

As a test of our code, in this section we investigate how 
well the slow-roll formulae work in the slow-roll regime for 
one of the usual potentials that do not predict large run- 
ning. As an example we will consider a simple quadratic 
potential, V = m^0^/2, which is a classic example of 
chaotic inflation. The second panel from the bottom in 
Fig. 121 effectively shows the dependence of z" /z on the 
number of the e-folds for inflation with such a potential. 
The behavior is monotonic and very smooth, which is 
due to the smoothness of the derivatives of the potential. 
Since z" /z scales as 2a^H^, we plot the quantity 



instead (compare to equation H12|l). 

The top two panels show the dependence of £//, r]H 
and on the number of e-folds. The only significantly 
non-zero term is en, which gradually grows to 1 at the 
end of inflation. The values of r]H and are typically 
smaller by roughly 10'^ and lO'' respectively. 

Figure |21 shows the primordial power spectrum pro- 
duced by the quadratic potential. The second panel from 
the bottom describes the error produced by the slow-roll 
approximations. The first order approximation gives less 
than 0.2% error in the observed range of wavenumbers k. 
The second order approximation works slightly better; 
the error is just above 0.1%. Both of these numbers are 
likely to be good enough for the upcoming experiments. 
This is because the accuracy at large scales is limited by 
the finite number of modes, while at small scales it is 
limited by the nonlinear evolution. So, while the overall 
amplitude could in principle be determined to an accu- 
racy of 0.1% when CMB and lensing information is com- 
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FIG. 2: Panels from the bottom to the top: 1. Potential 

V = m^<^^/2 in the range of cjfi's where the inflation occurs 
around 50 e-folds before the end of inflation. The scale in 

V is not COBE-normalized on this plot. 2. The dependence 
of z" /z on the number of the e-folds during the inflation. 
Number of e-folds A'' = corresponds to our arbitrarily chosen 
pivot scale of = 0.05/Mpc. 3. Plots of Hubble slow-roll 
parameters eH, fjH, ^h- The value of en is non-negligible 
whereas rjn and are essentially zeros. 4. Plots of horizon- 
flow slow-roll parameters ei, £2, £3 and the product of £263 for 
the same inflationary model. Since in this case en S> r]H,^H, 
the value of ei is essentially ei « en and £2 ~ £3 ~ 2eH. 
Nothing unexpected is going on here for this model. 



bined, it is unlikely that such a precision will be achieved 
separately at two widely separated length scales. 

Taking a more careful look at the error plot, one sees 
that the error curve in the observed area is basically a 
straight line, meaning that the main source of error is 
not the imprecise value of as but the error in Ug. Let 
us estimate now how precisely we need to know to 
get an error of, say 0.2%, in the observed range. The 
imprecision Sus in Ug will give us the uncertainty 
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0.002. 



(19) 



Taking the observed range of /c's to be from 10~^ Mpc~^ 
to 1 Mpc~^, we find that one needs to find Ug with the 
precision of Sus = 6-10^''. 

The same allowed uncertainty Sas in as is estimated 
from 
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—(5a, — In 
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= 0.002. 



(20) 



Therefore Sag — 3 • 10^^ is the error which we can make 
in determining as in order to get an error in the power 
spectrum of 0.2% at the edges of the observed range of 
/c's. 

The two top panels of Fig. |21 compare the numerically 
found dependence of and as on k to the one found 
from the slow-roll approximation with the first order as 
and either the first or second order for Us- One should 
compare the discrepancies between these to the values 
of 6ns and Sag. The characteristic value of is .964 
and the discrepancy between the exact value and the one 
found from the slow-roll approximation is comparable to 
Sus- Running as takes values around —6.5 • 10^"*. The 
discrepancy between the exact and the slow-roll values 
is very small in comparison to Sas- One can also notice 
that in this case |as| ~ 2Sas. Thus, even if we assigned 
as = 0, we would not get a significant error in the approx- 
imation of the primordial power spectrum of the scalar 
perturbations. 

To summarize this section, for standard inflationary 
potentials, the slow-roll approximation suffices even at 
first order when compared to the expected accuracy of 
existing and future experiments. The second order ap- 
proximation, while improving the accuracy, is not really 
necessary. The main error of slow-roll when considered 
in contrast to the numerical solutions is the inaccuracy 
in the slope n^; inaccuracies in higher order expansion 
terms, such as the running, are less important and can 
even be ignored. 



IV. POTENTIAL WITH A BUMP IN THE 
SECOND DERIVATIVE 



We want to construct a potential which will give us a 
strong running and crossing of the point rt^ = 1 in the 
observable power spectrum. We want to have > 1 at 
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FIG. 3: Panels from the bottom to the top: 1. Power spec- 
trum produced by inflation with potential V — m^(f>^/2. 
2. Errors made by approximating the true scalar power spec- 
trum by 1131 with calculating Us and Us from flrst order slow- 
roll formulae 1161 1171 , second order slow-roll formulae 1141 115II 
and calculating them through numerical derivatives. Numeri- 
cal third order takes into account third logarithmic derivative 
of the power spectrum /3s in parameterization 1131 . 3. The 
evolution of — 1 is calculated numerically and by the first 
and second order slow-roll formulae. 4. The evolution of Os 
is calculated numerically and by the slow-roll formulae (flrst 
and second order slow-roll are the same for Os). The error 
plot clearly shows that the main error comes from the impre- 
cision of Ws — 1, whereas the approximation for the as works 
well enough. One can also see that in the case of this poten- 
tial both first and second order slow-roll formulae for — 1 
overestimate the real value of rZs — 1. 



earlier times in inflation, while at later times we want to 
have Us < 1- To get the desired result, one can take two 
different potentials producing such features and smoothly 
connect them. 

One can rewrite slow-roll formulae H14I15() through the 
potential slow-roll parameters as 



ris — 1 = -6ey + 2t]v, 

as — WevW — 24ey — 2^v 

Now let us just choose our potential to be 

V{(l)) = 1-0.010-1.200^ 



(21) 
(22) 

(23) 



for all > 0. This choice provides about 50 e- folds of 
inflation after (p — 0. Since the local properties of the 
power spectrum are mostly determined by the local "his- 
tory" of the slow-roll parameters at the moment of the 
horizon crossing, we can get a red tilt of the scalar power 
spectrum w 0.80 in the area where the "history" be- 
fore point (p — is not very important. To get an ap- 
proximately symmetric shape of the power spectrum we 
choose V{(j)) to be 



V{(f)) = 1-0.010+1.2002 



(24) 



for all < 0. In this case for wave modes which cross the 
horizon far before the moment when the scalar field takes 
the value of = 0, the spectral index of the primordial 
power spectrum has a blue tilt ris ~ 1.20. Thus between 
these two regions the spectral index changes from 1.20 to 
0.80. We can unite formulae and 1(231) into 



F(0) = 1 - 0.010 - 1.2002 sign 0. 



(25) 



This potential has continuous first and second deriva- 
tives, but has a bump in its third derivative. This makes 



in H12|l discontinuous around = 0. According 

to [331 this produces oscillations in the power spectrum, 

which we can indeed see for the potential (|25|l . To avoid 

the oscillations we smooth out the sign function, chang- 
2 

ing it to — arctan(2OO0). In this case (|25|l changes to 



1/(0) = 1 - 0.010 - 1.2002- arctan(2OO0), (26) 

TT 

which is shown on Fig. ^ We have chosen 200 as the 
coefficient in front of in the arctan function so that the 
produced power spectrum has a nice shape as in Fig. [5| 
The two top panels of Fig. ^ show the behavior of the 
slow- roll parameters eu, rju, and ei, e2, ea, £263 cor- 
respondingly. While nothing unexpected happens to the 
behavior of the conventional Hubble slow-roll parame- 
ters en, TjH and there appears to be a singularity for 
the horizon-flow parameter £3. However, notice that the 
product e2e3 behaves smoothly and remains small due to 
the fact that the parameter £2 is changing its sign and 
therefore crossing through zero. Thus the parameteriza- 
tion of equation (jSJ introduces a singularity which is not 
physically present in the model. 
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FIG. 4: Panels from the bottom to the top: 1. Potential 1)26^ . 
2. Effectively this plot shows the dependence of z" j z on the 
number of the e-folds during the inflation. 3. Plots of Hubble 
slow-roll parameters e_f/ , r\H , £.h ■ Though the potential has a 
singular behavior, all the flow parameters are smooth. 4. Plots 
of horizon-flow slow-roll parameters ei, £2, £3 and the product 
of £263 for the same inflationary model. While everything is 
flne with £i, £2 and £2£3, the value of £3 indeed flips over 
inflnity. 



Figure shows the power spectrum produced by the 
model of inflation with the potential H25|) . The second 
panel from the bottom shows the errors made by different 
approximations. We again observe a similar picture for 
the slow-roll formulae. The main source of error for either 
the first or second order approximations comes not from 
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FIG. 5; Panels from the bottom to the top: 1. Power spec- 
trum produced by 12611 . 2. Errors made by approximating 
the true scalar power spectrum by (1131 with calculating Us 
and as from flrst order slow-roll formulae (llfcilllTII . second or- 
der slow-roll formulae 114II15II and calculating them through 
numerical differentiation. Numerical third order takes into 
account third logarithmic derivative of the power spectrum 
/3s in parameterization (I13II . 3. The evolution of — 1 is cal- 
culated numerically and by the first and second order slow-roll 
formulae. 4. The evolution of is calculated numerically and 
by the slow-roll formulae (first and second order slow-roll are 
the same for Us). The error plot clearly shows that the main 
error comes from the imprecision of — 1, whereas the ap- 
proximation for the Os works well enough. We see that in the 
case of this potential, flrst order slow-roll formula for — 1 
underestimates the real value, while the second order formula 
overestimates it. 
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the value of as but from the error in the value of n^. 
From the second panel from the top, we can estimate 
that the discrepancy is of the order of 0.01 for 
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which gives an error of 



1 h 
ons — in 

9 h ■ 



4% 



(27) 



(28) 



in the produced power spectrum at the edges of the 
observed range. Both the first and second order slow- 
roll approximations for Us work somewhat unsatisfactory. 
The first order slow roll underestimates Ug and the sec- 
ond order overestimates it by about the same amount. 

On the other hand, if our goal is to focus on running 
alone regardless of the slope and use just that property to 
deduce something about the potential, then the slow-roll 
does very well, since the differences between the slow-roll 
and numerical value of running are very small even at 
the lowest order in slow-roll. Extra terms in the expan- 
sion 1)13(1 further improve the accuracy. Adding running 
of the running improves the accuracy over the observed 
range from 1% to 0.2%. 

In summary, for potentials that lead to large running, 
slow-roll does not estimate the slope Ug very accurately 
at either first or second order, while the accuracy of the 
running as sufRces for the existing and future experi- 
ments. If we observe over a wide range of scales then it 
is useful to add the cubic term. Second order slow-roll 
does not seem to improve the accuracy. 



nonperturbed Hamilton- Jacobi inflationary equation 1^, 
which gives us an attractor solution i?(0) which in the 
inflationary class of problems almost does not depend on 
the initial condition. 

By following the method prescribed by one avoids 
solving the main attractor inflationary equation f^l, as 
pointed out by 38]. Indeed, the assumption ^S,h = for 
&\\i>M + l requires that H'^^\<j)) = for al\i>M + 2. 
Consequently, H{4>) is a polynomial of order M + 1: 

H{^) =Ho{l + Ai(b + A2<i? + A3<^3 + • • • + AM+i<t>^'^^). 

(32) 

In this case the function is an attractor solution of 

equation J^l with a potential in the form 



'pi 



327r2 



327r2 



Hi 



(v4i + ... + (M + l)^M+i</'''0 



- — (1 + ^1-; 



m: 



(33) 



Thus the only differential equation one needs to solve 
in order to match up the number of e-folds and the value 
of the scalar field 6 is 



dN 2^ 1 



dtp mpi yJKcj)) 
Here again e(0) is defined as in the equation ||SJ): 



(34) 



V. FLOW EQUATIONS SIMULATIONS 

Kinney (4Qj introduced a formalism based on the so- 
called flow equations, further discussed in |3l|. The ba- 
sic idea is that if one fixes the Hubble slow-roll parame- 
ters (jni at some point in time for e/f, rjH and ^S^h up to 
i = M and assumes that all the other Hubble slow-roll 
parameters are small enough that one can neglect them 
in one's calculations (i.e. ^£^h = for all £ > M -|- 1) 
then, without any other assumptions about inflation be- 
ing slow-roll, one can flnd the Hubble slow-roll param- 
eters at any other moment of time using the following 
hierarchy of linear ordinary differential equations: 



de 

dN 
drj 

dN 
d^ 
dN 



[^£-(£-1)77] 



(29) 
(30) 

(31) 



for aU £ = 2 . . . M assuming = 0. 

Usually when we set up an inflationary problem, we 
choose a potential V{(j)) and then reconstruct the form 
of the Hubble parameter during inflation using the main 



'"-pi 
47r 



Ai + 2A2 



M 



1 + Ai(j) + A2(, 



■■■ + A 



M+l' 



AM+1 



(35) 



We do not have to numerically solve the hierarchy of 
M differential flow equations. Instead we have analytical 
expressions for F(0) and H {(/>). 

The late attractor e = = and rj — const, found 
by 01' corresponds to the situation where the inflation 
proceeds to the value of the scalar field (p, which is a 
solution of the equation e — 



Ai + 2A2(f> + • • • + (M + l)AM+i(f> 



M 



0. 



(36) 



At this point if 7^ 0, then 77 ~ const 7^ due to the 
definition of f], which does not involve e at all. All the 
other = since any of them is a product of the first 
derivative of the Hubble parameter (which is zero) with 
some higher order derivatives. 

Peiris et al. [l^ made M — 9-th order flow equation 
simulations; about 40, 000 are shown as black dots on 
Fig. Fifty points fall into the range jn^ — 1| < 0.05 
and as < —0.02; these are shown in yellow. Among 
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FIG. 6: Panels from the bottom to the top: 1. Poten- 
tial reconstructed from one of the flow equations simulations. 

2. Effectively this plot shows the dependence of z" / z on the 
number of the e-folds during the inflation for this potential. 

3. Plots of Hubble slow-roll parameters en, rjH, (,h- All of 
them are smooth. 4. Plots of horizon-flow slow-roll parame- 
ters £1, £2, £3 and the product of £263 for the same inflationary 
model. While everything is fine with ei, £2 and £2£3, the value 
of £3 again flips over infinity. 



these point we have chosen 13 which fall into the narrow 
interval jris — 1| < 0.02, and we have reconstructed the 
corresponding inflationary potentials for the inflationary 
models which give such significant running, together with 
Us extremely close to 1. 

The bottom panel in Fig.^lshows a potential from such 
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FIG. 7; Panels from the bottom to the top: 1. Power spec- 
trum produced by potential on Fig. |S| 2. Errors made by 
approximating the true scalar power spectrum by (1131 with 
calculating and as from first order slow-roll formulae llfcil 
I17I I. second order slow-roll formulae 1141 1151 and calculat- 
ing them through numerical differentiation. Numerical third 
order takes into account third logarithmic derivative of the 
power spectrum l3s in parameterization (I13II . 3. The evolu- 
tion of Tis — 1 is calculated numerically and by the first and 
second order slow-roll formulae. 4. The evolution of Us is cal- 
culated numerically and by the slow-roll formulae (first and 
second order slow-roll are the same for as). The error plot 
clearly shows that the main error comes from the imprecision 
of ns — 1, whereas the approximation for the as works well 
enough. One can also see that in the case of this potential 
first order slow-roll formula for Us — l underestimates the real 
value, while the second order formula overestimates it. 
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a model with an unusually high value of the running as . 
We notice that there is a small dip in the potential. Some 
of the potentials with high as from the simulations had 
unrealistically high values of the tensor to the scalar ra- 
tio, but all of them had quite similar shapes. The second 
from the bottom panel of Fig. shows the characteris- 
tic behavior of the function z" / z which influences the 
scalar power spectrum as we have seen earlier. The top 
two panels show the dependence of the slow-roll param- 
eters on the number of e-folds. As in the other case with 
large running, we find a singularity for the horizon-flow 
slow- roll parameter £3, while the product £2^3 behaves 
smoothly and £2 crosses zero. 

The bottom panel in Fig. [3 shows the power spectrum 
of scalar and tensor perturbations produced by inflation 
with the potential under consideration. The second from 
the bottom panel shows the error produced by every one 
of the approximations for the power spectrum. We again 
see that both first and second order slow-roll formulae do 
not give a satisfactory result for Us- One of them again 
overestimates Hs] the other underestimates it. The error 
for either of the approximations is about 2-4%. 

The error introduced by the approximate formula for 
the running as is a bit smaller than the one for Us, but it 
is somewhat larger compared to the quadratic potential 
we considered in the previous section. 

Chen et al. 0| perform a similar analysis of slow-roll 
approximation. Using the flow-equations technique, they 
found discrepancy of larger than 0.01 for Us between sec- 
ond and third order slow-roll approximations for some of 
the models. Based on this fact they conclude that third 
order slow-roll is better. For the model we considered in 
this section we have found that the third order slow-roll 
does not improve the results of the second order approx- 
imation. In our calculations both formulas give identical 
results leading to approximately the same order of error 
as the first order approximation. 



VI. IS TRUNCATED TAYLOR EXPANSION 
GOOD? 
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FIG. 8: The bottom panel shows the numerically calculated 
power spectrum of the potential we used in section|3and it's 
approximations by the second order Taylor expansion and by 
Abazajian, Kadota and Stewart (AKS) parameterization 1381 1. 
The top panel shows the error produced by each of the pa- 
rameterizations. We have also added the error produced by 
the third order Taylor expansion here. AKS parameteriza- 
tion acceptably describes the true power spectrum in a very 
narrow range of fc's around the pivot point k — 0.05 Mpc~^. 
It does not improve the truncated Taylor expansion over the 
wider range of fc's. 



Recently Abazajian, Kadota and Stewart |42 
gued that if 



, ln{k/k^)\ > \ns 



II 



have ar- 



(37) 



then the traditional truncated Taylor series parameteri- 
zation is inconsistent, and hence it can lead to incorrect 
parameter estimations. One can notice that Taylor ex- 
pansions P{x) — ^fliZ* of functions or cos a; around 
a; = also violates the condition ai > a2X, but no one 
argues that these expansions are not valid. Abazajian et 
al. propose to use the parameterization 



InP(fc) =lnPo 



{ns - IT' 



Cts 



(38) 



instead. 



There is one significant disadvantage of this approach. 
In particular, using this parameterization to describe V 
as a function of k, one is able to describe only a growing 
or decreasing function, which can be easily seen from 
the form of the function. The models we study in this 
paper produce scalar power spectra which are not purely 
growing or decreasing (e.g. see Fig.O. 

In the previous section we have considered the po- 
tential which produces power spectrum satisfying equa- 
tion H37() . On Fig. 121 we compare the traditional trun- 
cated to second and third order Taylor expansion and 
the parameterization H38|l . We find that the parameteri- 
zation H38|) gives a significantly larger error than, e.g. the 
second order Taylor expansion. 

Thus we found that in this particular case though equa- 
tion (|37|l holds, truncated Taylor expansion is a good ap- 
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proximation and the AKS approach does not improve it. 
There might be models for which equation (|38|l works 
better than Taylor expansion, but it is definitely not an 
improvement for a general case and should be used with 
caution, if at all. 

VII. CONCLUSIONS 

In this paper we have explored the accuracy of the 
slow-roll approximation given the observational con- 
straints on the primordial scalar and tensor power spec- 
tra. The current constraints can be roughly described 
by the tensor to scalar ratio r < 1, small deviation 
from the scale invariance of the scalar power spectrum, 
\ns — 1| < 0.05 and small but possibly nontrivial running, 
|as| < 0.03. These constraints allow for the particular 
case where ~ 1 and ja^j > 0.01, which has previ- 
ously been argued to not satisfy the slow-roll condition. 
We have computed exact numerical solutions for the con- 
sidered potentials and compared them to those obtained 
from the first and second order slow-roll approximations. 

We have found that for the potentials explored here, 
there is no substantial difference when using first or sec- 
ond order slow-roll formulae for the power spectrum in- 
dex Us . Both of them either work well in the case of small 
running or have a comparable error in the case of non- 
negligible running. Adding extra (cubic in In fc) terms in 
the approximation for the scalar power spectrum extends 
the accuracy to a larger range of scales, but this accuracy 
is most likely not necessary for existing and near future 



experiments. If the values of and are known with 
the precision Sug — 6 • 10^^ and Sag — 3 • 10^^, then the 
scalar power spectrum will have an error of about 0.2% 
at the edge of the observable range of wavenumbers k's. 

The horizon- flow basis e„_|_i = dln\en\/dN introduces 
an artificial singularity for inflationary models with neg- 
ative running and the value of the spectral index crossing 
1. Such a divergence in one of the horizon- flow param- 
eters does not indicate that the slow-roll approximation 
has been badly broken. We find that the slow-roll is still 
accurate at the 1-2% level and most of the error comes 
from inaccuracies in the evaluation of the slope itself, and 
not the running. Thus the first order slow-roll approx- 
imation is sufficiently accurate for the current observa- 
tions. Only if the running turns out to be large, while 
the slope remains close to scale-invariant, are exact nu- 
merical calculations required to achieve sub-percent ac- 
curacy. In the appendix we present a short guideline on 
performing such calculations. One can request the code 
directly from the author. 
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APPENDIX: INFLATIONARY EQUATIONS 

In this appendix we describe the technical details of 
the code we ran to get the results presented in the main 
part of the paper. The code is given a potential V{(j)) 
and some point (jjQ which lies in the observable range of 
wave-modes and, say, corresponds to the moment when 
wavelengths with k = 0.05 Mpc~^ exit the horizon. We 
want to find the power spectrum produced by inflation 
with the potential Vi^cj)). For this purpose we first have 
to go backwards in time about 50 e-folds and then start 
the infiation there. This guarantees that the infiation- 
ary dynamics are not affected by the choice of the initial 
condition and we indeed have the attractor solution. 

Now we evolve the universe from our "beginning of in- 
flation" to the end of inflation, the moment which is de- 
termined by the violation of the inequality a > 0. This 
part is described below in the "non-perturbed inflation- 
ary equations" section. Usually we require 50 to 70 e- 
folds between (/jq and the end of inflation. 

After we already have the complete background history 
of the evolution of the universe during the inflationary 
stage of the expansion, we can start working out the evo- 
lution of the perturbations during inflation, as discussed 
in the second part of the appendix. 



1. Non-perturbed inflationary equations 

The unperturbed dynamics of inflation are described 
by the equation of motion of the scalar field (j) with po- 
tential V{(j>) in the expanding universe with the Hubble 
parameter H = d/a 



4> + 3H,j) + V'{<l)) =0 



(A.l) 



and the Friedman equation with only the scalar field com- 
ponent present in the universe 



Stt 
3^ 



(A.2) 



The equations (|A.1I IA.2|) are equivalent to the pair of 
Hamilton- Jacobi equation ^ and 



47r 



(A.3) 



The Hamilton- Jacobi equation connects the Hubble 
parameter and the value of the potential of the scalar 
field during the inflation. In the case when we know 
the behavior of the Hubble parameter it is easy to find 
the potential. The method of flow equations is entirely 
based on this fact. In contrast, if we know the shape of 
the potential and want to reconstruct the behavior of the 
Hubble parameter, the problem is not as simple. First of 
all, as for any first order differential equation, we would 
like to have an initial condition Hq = iJ(0o). Due to 
the attractor nature of the equation its solution does 
not really depend on the initial condition Hq (we have 
found from numerical simulations that one needs about 
6 e-folds to forget the history). Thus it does not really 
matter which initial condition we choose. 

Hamilton-Jacobi equation requires that 



(A.4) 



If we are going to use a method such as Runge-Kutta for 
the integration of the differential equationlj^l, we might 
try values of H which would violate the inequality (jA.4|l . 

To avoid this complication, we reparametrize our equa- 
tion using a new function (5(0) so that 



(A.5) 



Then substituting our new definition into equation 
we get 



H'{ct>) 



,5(0)/2 



(A.6) 



Combining this with the expression for H' obtained from 
the direct differentiation of H in (|A.5p . we get a differ- 
ential equation for (5'(0) 



S' = -Vl 



V 



TOpl 



(A.7) 
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This equation is much more pleasant to deal with numer- 
ically than equation (Q, since it does not have a weird 
boundary for 5, as H did before. One can also check the 
attractor nature of the equation (|A.7p . that it does not 
remember the prior history. We see now that in the case 
when the potential is changing slowly S' and we have 
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ml,iV'/Vf 
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We can use this approximate solution of the equation 
as the initial condition for our differential equation since 
it is quite close to the true solution and it will make 
our numerical solution evolve into the attractor solution 
faster. 

One can check that the expressions for en, ijh and ^^h 
are given by the following formulae 



m: 



An \H 

2 TT// 



m 



Att H 



3 V 



V ^e«(l + e^)' 



(A.IO) 
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3(e + 7ii) — rf 



3m2j V" 
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From these expressions we can expect that in general e 
and rj are continuous functions, whereas does not have 
to be continuous at points where V" is not continuous. 

The condition for inflation to take place (a > 0) follows 
from the derivative of the Friedman equation 
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(A.12) 
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V{(j3){l-2e^). (A.13) 



The first of these two equations implies that the end of 
inflation happens when the inequality e < 1 is violated. 
The same thing occurs when the inequality (5 < — In 2 is 
violated in the second equation. The latter also means 
that the inflation continues while the kinetic energy of 
the inflaton is less than half of its potential energy 



K 

n 



1/2 



J{4>) < }_ 



(A.14) 



Thus we come to a physical definition of our parame- 
ter e'^^''^^ as the ratio of the kinetic energy 0^/2 to the 
potential energy V{4)). 

In the next subsection we will be working with infla- 
tionary perturbations and it will not be very convenient 
for us to work with the value of the scalar field (j) as an 
independent variable. For this purpose we will use the 
number of e-folds defined as 



N = In 



{aH)_ 



(A.15) 



Note that this is the actual number of e-folds and is not 
the same as iV = ln(a/ao). The connection between TV 
and (f) is determined through the derivative 



dN 20F 1 - e(0) 



(A.16) 



rf0 rupi y/e{(f)) 
We are almost done describing the background evolu- 



tion of the universe, except we have not yet chosen the 
initial value of the scalar field We only have the value 
(A. 9) '/'o which corresponds to the moment when the mode 
k — 0.05 Mpc~^ exits the horizon. We want to move 
backwards in time for about 50 e-folds. Equation l|A.7|l 
has an attractor behavior only when we are moving in the 
positive direction along the (/>-axis. It diverges from the 
attractor solution in the negative direction. As a useful 
trick, let us modify equation |^ to the following form 
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V{c^)~H\<i>) 



(A.17) 



(A.ll) 

In this form, when we move backwards in time the value 



of is bound by the value of ^S-KV/Zm^^ from the 

top and the solution cannot diverge. In addition we tem- 
porarily redefine (5(0) to satisfy 



S^pi 



(A.18) 



Thus we get an equation analogous to the equation ljA.7|) 



(5' = vT^ 



V 



4V37r 



(A.19) 



Equation HA.19|I does not carry any physical mean- 
ing; we just use this equation to go "upwards" to the 
higher values of the potential, still tracking the general 
behavior of V{4)). If we go backwards in time 50 e- 
folds using HA.19|I and then forward in time 50 e-folds 
using (|A.7p . we will not return to the same point 
since the behavior of 6{(t>) in the equation (|A.19p is de- 
termined by the area which is to the right of the current 
value of (f) and in the equation (|A.7|I is determined by the 
area which is on the left side. Nevertheless, this method 
gives us a good estimate of what initial value of (j)i we 
should take. 

It is also worth mentioning that this approach is not 
more difficult to deal with than the inflationary flow 
equations (|29H31|1 . 
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2. Perturbation equations 

a. Scalar mode 

The algorithm for finding the scalar mode primordial 
power spectrmn is described in the main text (see equa- 
tion H10|l and below). Here we will just mention some 
technical details. 

Equation (|l()|l is not very convenient to solve in its 
current form. First of all we would like to set the inde- 
pendent variable, the conformal time r, in such a way 
that T ^ — as inflation goes on. In this case we would 
be able to numerically integrate equation (|10|l up to as 
small values of t as we want. But in numerical realiza- 
tions we cannot really choose such an initial value of Ti 
that gives us r — > —0 at the end of the inflation. Suppose 
that at the end of the inflation we have r — > 1 — 0. In 
this case the numerical error on t will be of the order 
of 10~^^ which is a reasonable machine precision. Hence 
the limit on corresponding dr is of the same order and 
we can explore the range of changing the scale factor a 
from ^ 1 to ~ 10^*^, i.e. about 35 e- folds. This might be 
enough, but to be safe we will use a different indepen- 
dent variable, the true number of e-folds N defined by 
equation IjA.lSjl which is the same as 
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Then the mode equation 1)1 Q|l can be rewritten as 
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where coefficient a, b and c can be exactly expressed 
through en, rjn and ^S,h as 



b = -2e- 



c = e 77 

2 ' 



+ 2677, 

- 2eri 



2^^ 



2 ^ 



(A.22) 
(A.23) 

(A.24) 



In equation l|A.21l) . k is the wavelength of the interest, 
while fco and A'o are constants conveniently chosen for 
normalization purposes. 

Further, equation 11U|) has a solution 
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at the beginning of the inflation when t 
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k^ ^ — —. We also know the approximate behavior of 
z ar^ 



Uk at later times when r —0 and fc^ ^ 
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Thus it is natural to decompose Uk into growing and 
oscillating parts 
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where both functions A and 4> are real functions of con- 
formal time r or of the true number of e-folds N. Then 
the equation l|A.2ip can be split into 4 ordinary differen- 
tial equations with 2 new functions Ap and (j)p defined as 
below 
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This system of differential equations looks a bit more 
complicated than the single equation (|10|) , but it is actu- 
ally much easier to solve numerically. Indeed, at earlier 
times we have dA/dN = Ap — 0. This instaneously gives 
us the initial condition on d(i>/dN = (f)p from (|A.29|I 
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as T — > — oo, i.e. N — > —oo. To be consistent with the 
initial condition on 
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as T ^ — oo we also require that 
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-oo. As the inflation continues, the terms 
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and 0p will balance each other on the right hand side of 
the equation (|A.29p until Ap is not negligible in compar- 
ison to 1 in equation (|A.3ip . Thus, around N = ln(A:/A:o) 
the oscillating part 4>p will decrease more rapidly than be- 
fore, finally exponentially dropping to zero. At the same 
time Ap, and therefore A, start exponentially growing. 
The final power spectrum is given by 
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Thus we even do not need information about the phase 
(f) and we can freely drop equation (|A.30p from our sys- 
tem. Also while being in the stage of inflation where Uk 



14 



has an oscillatory behavior, if one were to use the usual 
method without our substitution, one would have to find 
the values of for at least 6 points per oscillation pe- 
riod. However with our substitution, we easily pass this 
area, which does not have any interest for us since we 
analytically know the behavior of Uk here, and therefore 
move directly to the place where we cannot solve it ana- 
lytically. By our estimates this technique gives a gain of 
a factor of 10 in computational time, which is of particu- 
lar interest if one wants to calculate the power spectrum 
for e.g. 100 wavemodes. 



with the same initial condition 

as r — !• — oo, where a is the usual scale factor of the 
Friedman universe. One can show that 



h. Tensor mode 

The calculation of the tensor mode power spectrum 
of perturbations is absolutely analogous to the one for 
scalars, except instead of equation H10|l one has to solve 



Mode equations for the amplitude and the phase of the 
wave (jA.28IA.3ip of the tensor mode look similar except 
in the equations (|A.22IA.24p where we have to change c 
to d defined as 



1 



-€. 



(A.36) 



